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Abstract 

We have tested the leading correction-to-scahng exponent w in 
0(n)-syninietric models on a three-dimensional lattice by analysing the recent 
Monte Carlo (MC) data. We have found that the effective critical exponent, 
estimated at finite sizes of the system L and L/2, decreases remarkably within 
the range of the simulated L values. This shows the incorrectness of some claims 
that uj has a very accurate value 0.845(10) at n = 1. A selfconsistent infinite 
volume extrapolation yields row estimates w « 0.547, uj « 0.573, and uj « 0.625 
at n = 1, 2, and 3, respectively, in approximate agreement with the correspond- 
ing exact values 1/2, 5/9, and 3/5 predicted by our recently developed GFD 
(grouping of Feynman diagrams) theory. We have fitted the MC data for the 
susceptibility of 3D Ising model at criticality showing that the effective critical 
exponent r/ tends to increase well above the usually accepted values around 
0.036. We have fitted the data within [I/;8L], including several terms in the 
asymptotic expansion with fixed exponents, to obtain the effective amplitudes 
depending on L. This method clearly demonstrates that the critical exponents 
of GFD theory are correct (the amplitudes converge to certain asymptotic val- 
ues at L — > c»), whereas those of the perturbative renormalization group (RG) 
theory are incorrect (the amplitudes diverge). A modification of the standard 
Ising model by introducing suitable "improved" action (Hamiltonian) does not 
solve the problem in favour of the perturbative RG theory. 

Keywords: Xip'^ model, Ising model. Binder cumulant, Monte Carlo data, crit- 
ical exponents 

1 Introduction 



Since the exact solution of two-dimensional Ising model has been found by On- 
sager [|l|, a study of various phase transition models is of permanent interest. Nowa- 
days, phase transitions and critical phenomena is one of the most widely investigated 
fields of physics. Remarkable progress has been reached in exact solution of two- 
dimensional models [Q]. Recently, we have proposed |Q] a novel method based on 
grouping of Feynman diagrams (GFD) in ip'^ model. Our GFD theory allows to 
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analyse the asymptotic solution for the two-point correlation function at and near 
criticality, not cutting the perturbation series. As a result the possible values of ex- 
act critical exponents have been obtained ||3| for the Ginzburg-Landau (ip^) model 
with 0(n) symmetry, where n = 1, 2, 3, ... is the dimensionality of the order param- 
eter. Our predictions completely (exactly) agree with the known exact and rigorous 
results in two dimensions Q, and are equally valid also in three dimensions. In 
we have compared our results to some Monte Carlo (MC) simulations and exper- 
iments 0, |5|, ^. An additional comparison to MC data has been made in P]. It 
has been shown ^ ^ that the actually discussed MC data for 3D Ising @, and 
XY [S] models are fully consistent with our theoretical predictions, but not with 



those of the perturbative renormalization group (RG) theory 11, 12|. Some data 
for 3D Heisenberg model also have been discussed in [^. However, these data, 
likely, are not accurate enough and here we reconsider the estimation of the critical 
point based on more recent MC results. From the theoretical (mathematical) point 
of view, the invalidity of the conventional RG expansions has been demonstrated 
in 1^]. The current paper, dealing with numerical analysis of the three-dimensional 
Xip^ and Ising models, presents one more confirmation that the correct values of 
critical exponents are those predicted by the GFD theory. 



2 Xcp model and its crossover to Ising model 

Here we discuss a ip'^ model on a three-dimensional cubic lattice. The Hamiltonian 
of this model, further called Ac^^ model, is given by 



H/T = J2 \ -2«^ J2 'P^'P^+f^ + V'x + ^ ("y 



V^x - 1 j I , (1) 

where the summation runs over all lattice sites, T is the temperature, 
(^x — oo;+oo[ is the scalar order parameter at the site with coordinate x, ft 
is a unit vector in the //-th direction, k and A are coupling constants. Obviously, 
the standard 3D Ising model is recovered in the limit A — > cxd where ip'^ fluctuations 
are suppressed so that, for a relevant configuration, — 1 oi" '/'x — il holds. The 
MC data for the Binder cumulant in this Xip^ model have been interpreted in ac- 
cordance with the e-expansion and a perfect agreement with the conventional RG 



values of critical exponents has been reported in [13|. According to the definition 
in [13 1, the Binder cumulant U is given by 



where m = L ^ X^x V'x is the magnetization and L is the linear size of the system. 



Based on the e-expansion, it has been suggested in |13| that, in the thermodynamic 
limit L —5- cxD, the value of the Binder cumulant at the critical point k = Kc(A) 
and, equally, at a fixed ratio Za/Zp = 0.5425 (the precise value is not important) of 
partition functions with periodic and antiperiodic boundary conditions is a universal 
constant U* independent on A. We suppose that the latter statement is true, but 
not due to the e-expansion. It is a consequence of some general argument of the 
RG theory: on the one hand, U is invariant under the RG transformation and. 
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on the other hand, an unique fixed point (not necessarily the Wilson-Fisher fixed 
point) exists in the case of an infinite system, so that U = U* holds at L — > oo and 
K = Kc{\) where U* is the fixed-point value of U. The above conclusion remains true 
if we allow that the fixed point is defined not uniquely in the sense that it contains 
some irrelevant degree(s) of freedom (like c* and A in the perturbative RG theory 
discussed in Sec. 2 of Q) not changing U. The numerical results in confirm the 
idea that liniL^oo U (L) = U* holds at criticality, where U* is a universal constant 
independent on the specific microscopic structure of the Hamiltonian. 

3 Estimation of the correction exponent lu 

Based on the idea that U* is constant for a given universality class, here we estimate 
the correction-to-scaling exponent u>. According to [0, corrections to finite-size 
scaling for the magnetization of the actual 3D Ising and Xip'^ models are represented 
by an expansion in terms of L~'^ where uj = 1/2. One expects that the magnetization 
(Binder) cumulant (§) has the same singular structure. Since lim^^oo U{L, A) = U* 
holds at a fixed ratio Za/Zp, a suitable ansatz for estimation of u> is |l3| 

[/(L,Ai) -f/(L,A2) -const •L^'^ at Z^/Zp = 0.5425 , (3) 

which is valid for any two different nonzero values Ai and A2 of the coupling constant 
A. The data for AU{L) = U{L, 0.8) - U{L, 1.5) can be read from Fig. 1 in Q (after 
a proper magnification) without an essential loss of the numerical accuracy, i. e., 
within the shown error bars. Doing so, we have evaluated the effective exponent 

cOeffiL) = In [AU{L/2)/AU{L)] / In 2 , (4) 

i. e., We//(12) ~ 0.899, We//(16) ~ 0.855, and We//(24) ~ 0.775. These values 
are shown in Fig. || by crosses. Such an estimation, however, can be remarkably 
infiuenced by the random scattering of the simulated data points, particularly, at 
larger sizes where AU (L) becomes small. This effect can be diminished if the values 
of AU (L) are read from a suitable smoothened curve. A comparison to the original 
results (without the smoothening) provides some objective criterion of the accuracy 
of such estimations. We have found that AU{L) within L G [7; 24] can be well 
approximated by a second-order polinomial in as shown in Fig. Without 

any claims about validity of such an approximation well outside of this interval, we 
can consider the least-squres fit in Fig. |^ as an appropriate smoothened curve from 
which we read We//(16) ~ 0.8573 and a;e//(24) ~ 0.7956. These values are depicted 
in Fig. |l] by empty circles. As we see, the results are similar to those obtained by 
a direct calculation from the original data points (crosses). However, the values 
obtained from the smoothened curve (circles) are more accurate and reliable. As 
regards the smallest size, we suppose that the original estimate ti;e//(12) — 0.899 
is accurate enough even without any smoothening, since the values of A[/(6) and 
AU{\2) are large relative to the statistical errors. 

In such a way, we see from Fig. || that the effective exponent u>f,ff{L) decreases 
remarkably with increasing of L. According to GFD theory, iVeff{L) is a linear 
function of L^^l'^ at L ^ cxo, as consistent with the expansion in terms of L^^ 
where uj = 0.5. More data points, including larger sizes L, are necessary for a 
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Figure 1: Effective correction-to-scaling exponent uJeff{L) in the 0(n)-synimetric Ai^'* 
model with n = I (empty circles and crosses) and n = 3 (triangles), and 0(2)-symmetric 
dd — XY model (solid circles) depending on the system size L. The linear least-squares 
fits give row estimates of the asymptotic oj values 0.547, 0.573, and 0.625 at n = f , 2, and 
3, respectively. The corresponding theoretical values of the GFD theory 1/2, 5/9, and 3/5 
(used in the L^"^ scale of the horizontal axis) are indicated by arrows. The dot-dot-dashed 
line shows the value 0.845(10) proposed in for the 31? Ising universality class {n = 1). 
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reliable estimation of the asymptotic exponent u = lim^^oo '^e//(^)- Nevertheless, 
already a row linear extrapolation in the scale of L~^/^ with the existing data points 
yields the result u ~ 0.547 which is reasonably close to the exact value 0.5 (horizontal 
dashed line in Fig. |l]) found within the GFD theory. The corresponding least-squares 
fit with circles (at L = 24, 16) and cross (at L = 12) is shown in Fig. |^ by a straight 
solid line. It is evident from Fig. ffl that the final result uj = 0.845(10) (horizontal 
dot-dot-dashed line) reported in represents some average effective exponent for 
the interval L € [6; 24]. It has been claimed in ||l^ that the estimates for lo (cf. 
Tab. 2 in |l^) are rather stable with respect to Lmin, where Lmin is the minimal 
lattice size used in the fit. Unfortunately, the analysis has been made in an obscure 
fashion, i. e., giving no original data, so that we cannot check the correctness of this 
claim. Besides, the estimates in Tab. 2 of [^] has been made by using an ansatz 

U{L, X) = U* + ci{X)L-^ at Za/Zp = 0.5425 , (5) 

which is worse than Namely, (P) and (|5|) are approximations of the same order, 
but (^) contains an additional parameter U* which is not known precisely. The re- 
sults of an analysis with the ansatz reflected in Tab. 5 of [13 1, are not convincing, 
since only very small values of Lmin (up to Lmin = 6) have been considered. 

In any case, we prefer to rely on that information we can check, and it shows that 
the claim in |l^ that io = 0.845(10) holds with ±0.01 accuracy cannot be correct, 
since LJeff{L) is varied in the first decimal place. 

We have made a similar estimation of uj for 0(n)-symmetric spin models, namely, 
for the dynamically diluted XY {dd — XY) model simulated in |jl^ {n = 2) and for 
0(3)-symmetric X^p'^ model simulated in [|l5[. In the case of the dd — XY model, 
parameter D (cf. Eq.(6) in |14]) plays the role of A in (p|). The data for the Binder 



cumulant in Fig. 1 of [|1^ look rather accurate, i. e., not scattered. This enables us 
to estimate Wg// just from the data at D = 1.03 and D = oo (XY model). The 
resulting values of i^e// are depicted in Fig. Q by solid circles. The scale of L~'^ is 
used, where w = 5/9 is our theoretical value of the correction-to-scaling exponent 
at n = 2 consistent with the general hypothesis proposed in 0]. As we see, the solid 
circles can be well located on a smooth line which, however, is remarkably curved at 
smaller sizes. Due to the latter reason, we have used only the last three points (the 
largest sizes) for the linear fit (solid line) resulting in an estimate ui ~ 0.573 which 
comes close to our theoretical value a; = 5/9 = 0.555 . . . 

The estimates of u^eff ^or the 0(3)-symmetric Xf'^ model are depicted in Fig. |l| 
(in the scale of L"^ with our oj value 0.6) by triangles. The data have been extracted 



from Fig. 1 of [15| at Ai = 2 and X2 = oo. The obtained oJeff{L) values at L = 12, 16, 
and 24 well lie on a straight line (tiny-dashed line in Fig. |^) yielding an asymptotic 
estimate a; ^ 0.625 which is reasonably close to our theoretical prediction w = 0.6. 
Hence, the ci;e//(32) value deviates upwards unexpectedly. We suppose, this is due 
to an inaccurate simulation of the largest size L = 32, as explained below. We have 
depicted in Fig. | the ratio R{L) = {U{L, 2) - U*)/{U* - U{L, 00)) evaluated from 



the data in Fig. 1 of [15| with U* = 1.14022 (the average over three estimates at a 
fixed Za/Zp given in Tab. 2 in |l^). According to (|5|), R{L) tends to a constant at 
L 00. The i?(24) data point well lie on the smooth line (parabola) representing 
the least-squares fit to four smallest sizes L = 6, 8, 12, and 16, whereas the simulated 
R{32) value drops down unreasonably. So, it looks like a wrong simulation has been 
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Figure 3: The ratio R{L) = {U{L, 2) - U*)/{U* - U{L, oo)) for 0(3)-symmetric model 
estimated from the MC data of (Tsj. The scale of with ui — 0.6 is used. The solid line 
(parabola) represents the least-squares fit including only four smallest sizes. 



made at L = 32 to confirm the known RG estimate uj ~ 0.8. Thus, our extrapolation 
in Fig. m (tiny-dashed line), omiting the point with L = 32, is justified. 

In summary, the extrapolated co values (Fig. |^) in all three cases n = 1,2,3 are 
reasonably close to our theoretical values 1/2, 5/9, and 3/5 indicated by arrows. 
Only a small systematic deviation is observed. This, likely, is due to the error of 
linear extrapolation: the uJeff{L) plots have a tendency to curve down slightly. The 
conventional (RG) estimate uj ~ 0.8 more or less corresponds to effective exponents 
for currently simulated finite system sizes, but not to the asymptotic exponents. 

The data for n = 4 also are available in |15|. Unfortunately, they are too much 
scattered for the actual analysis. 



4 The critical coupling of 3D Ising 
and Heisenberg models 

A conventional method to determine the critical exponent r/ is a fit of the suscep- 
tibility data at criticality. For this, however, we need a very accurate value of the 
critical coupling (temperature). In this section we discuss the estimation of the 
critical coupling Pc for 3D Ising and Heisenberg models. 



The critical point of the standard 3D Ising model has been estimated in |16] with 
a 7-digit accuracy, i. e. /3c = 0.2216545. We have made our own fits with the MC 



data of [16 1 to check the accuracy of this estimation, and have obtained the same 
value within error bars of 10"''. We will use in our further analysis also a similar 
estimate /3c = 0.383245 [16| for the so called "improved" 3D Ising model. 

The critical coupling of the classical 3D Heisenberg model is known much less 
accurately than that of the 3D Ising model. Some of the known estimates are 
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Figure 4: The Binder cumulant U vs system size L in the 3D Heisenberg model at = 
0.6929 (bottom), /3 = 0.692955 (middle), and fi = 0.693 (top). Symbols depict the MC 
data of whereas lines represent the least-squares approximations of these data by a 
parabola. The solid line, coinciding with the universal critical value of the Binder cumulant 
U* = 0.61993(3), corresponds to the critical coupling (3^ 0.692955. 



f3c = 0.6929(1) IT3], f5c = 0.693035(37) |8|, and fic = 0.693001(10) ]T9|. 

In principle, the location of the critical point can be found with a high accuracy 
and reliabihty by simulations of the Binder cumulant U (/3, L) in close vicinity of 
the critical point, as it has been done in [^] for the 3D Ising model. Taking into 
account the leading and the subleading corrections to scaling, we have 

C/(/3c, L) ~ [/* + aiL-'^ + a2L-'^'^ . (6) 

If the universal value of U* is known with a high precision, then the critical value 
of f3 at which U {(3, L) coincides with (^) is well defined. Namely, at /? = /3c the 
quadratic (least-squares) extrapolation of f7(/3, L) in the scale of L""^ should yield 
U* at —5- 0. In Fig. ^, we have shown the results of such an extrapolation with 
uj = 0.6 (our theoretical value) at three different values of i. e., 0.6929 (lower 
dashed line), 0.692955 (solid line), and 0.693 (upper dashed line). The data for 
U{(3, L) have been extracted from Fig. 1 in ||20|] via an approximation f7(/3, L) ~ 
uIPo,L) + 0.06L^-^(/3 - Po), where /?o = 1/1.4432. Note that, in distinction to (g), 
now we use the conventional definition U = 1— (l/3)(m^) / [m?)'^. A suitable estimate 



of [/*, taken from Tab. 2 in |15|, is then U* = 0.61993(3). This value with the error 
bars is indicated in Fig. ^ by an arrow. As we see from Fig. ^, the estimate /?c — 
0.692955 is consistent with this value of U* . Our estimation is rather stable, i. e., 
if the extrapolation is made in the scale of L~^'^, we get practically the same result 
— 0.692957. Taking into account the curvature of the U{f3, L) plot at /3 = 0.693, 
it is unlikely that 0.693 could be the correct value of (3c yielding U* = 0.61993(3) at 
L — > CX3. According to the strong variation of the extrapolated U* value with /?, it 
is plausible that the error of our estimation f3c ~ 0.692955 is about 0.00001 or even 
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smaller. 

Our result agree within the error bars with the value 0.6929(1) of while the 
error bars of the estimates = 0.693035(37) |T|1 and Pc = 0.693001(10) ll|, in 
our opinion, are underestimated. As regards the value of [18j, this is a result of a 
linear extrapolation (in the scale of 2^-i/0-7036-j q£ temperature values corresponding 
to extrema points for several physical quantities. However, the shift in /3 as large 
as 0.0001 is practically invisible in the scale of Fig. 3 in 1 18]. If one allows that 
the lines in this figure are curved very slightly, then even larger deviation in the 
extrapolated /?c value is possible. In other words, the proposed error bars ±0.000037, 
obviously, include only the statistical error and are underestimated since they ignore 
the possible systematic deviation due to the neglected corrections to scaling. A 
similar problem with the subleading correction persists in |1£]. Namely, a linear 
extrapolation in Fig. ^ would give a larger Pc value, close to 0.693 (in agreement 
with that of |0]), but this value is shifted down to 0.692955 due to the subleading 
correction in the asymptotic expansion of U. This correction has been neglected in 
the analysis of the Binder cumulant crossings in [19|. We have taken into account 
both the leading and the subleading corrections to scaling and, therefore, our value 
0.692955 is more accurate than those proposed in [18| and [19|, unless the data of |2C] 
contain large systematical errors. 



5 Fitting the susceptibility data at criticality 

In this section we discuss some fits of MC data at criticality. According to the 
finite-size scaling theory, the susceptibility x near the critical point is represented 
by an expansion 

X = L"--^ i^o{L/i) + E L--^gi{L/i)^ , (7) 

where gi{L/S^) are the scaling functions, ^ is the correlation length of an infinite 
system, rj is the critical exponent related to the k~'^~^'^ divergence of the correlation 
function in the wave vector space at criticality, and lvi are correction~to~scaling 
exponents, loi = uj being the leading correction exponent. The correlation length 
diverges like ^ oc t~'^ at t — > 0, where t = 1 — /3//3c is the reduced temperature. Thus, 
for large L, in close vicinity of the critical point where tL^^'^ ^ 1 holds Eq. (|^ can 
be written as 

X = aL^-^ (^l + J^biL-'^' +Sit,L)^ , (8) 

where a = go{0) and bi = gi{0)/go{0) are the amplitudes, and 5{t,L) is a correction 
term which takes into account the deviation from criticality. In the first approxima- 
tion it reads 

6{t,L) c^c-tL^/" , (9) 

where c is a constant. 

We start our analysis with the standard 3D Ising model with the Hamiltonian 

H/T = -PY,<^^<TJ . (10) 
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Figure 5: The fits of In (x/L^) data at criticality (ansatz (||)) shifted by a constant c. Circles 
represent the MC data for 3D Ising model jl^ at f3 — 0.2216545 (c — 0, empty symbols) and 
3D Heisenberg model at /3 = 0.692955 (c = 0.55, solid symbols). The corresponding fits 
with our (GFD) exponents (In a = 1.065289, 6i = -2.72056,62 = 8.18636,63 = -10.49614 
and Ina = 0.207324, 6i = —1.22546,62 = 1.85823) are shown by solid lines, whereas those 
with the exponents of jl|, |2l| (Ina = 0.430933, 6i = 0.05850,62 = -7.74767,63 = 12.42890 
and Ina = -0.150242, 6i = 0.03947,62 = -0.45033)_ by tiny-dashed lines. The empty 
boxes are MC data for 3-component 3D XY model [p2|, shifted by c = 0.85. 



The critical point of this model with a 7-digit accuracy is Pc — 0.2216545 (Sec. Q). 
From the maximal values of the derivative dln{'m?) /dfi = dlnx/d(3 evaluated in |23] 
we conclude that the shift of /3 by 10~^ produces the variation of Inx at L = 96 
near f3 = j3c, which does not exceed 4.7 • 10~^ in magnitude. The latter means that, 
with a good enough accuracy, we may assume that /3c is just 0.2216545 when fitting 
the susceptibility data at criticality within L G [4; 128]. Here we mean the MC data 
given in Tab. 25 of |l^. We have made and compared several fits of these data to 
ansatz (^) with 6{t,L) = (more precisely, to the corresponding formula for In^) 
for two different sets of the critical exponents, i. e., our (GFD) and that proposed 
in 1 13]. The fits made with our exponents systematically improve relative to those 
made with the exponents of [13|, as the system sizes grow and the approximation 
order increases. The necessity to include several correction terms is dictated by the 
fact that corrections to scaling are rather strong. According to the least-squares 
criterion, the fit with our exponents rj = 1/8 and lji = 1/2 becomes better than 
that provided by the more conventional exponents rj = 0.0358(4), loi = 0.845(1), 
UJ2 = 2oji, and uj^ = 2 |jl3[ starting with Lmin = 28 (i. e., L G [Lmin] 128]), if two 
correction terms (/ = 1,2) are included. In the case of three correction terms it 
occurs already at Lmin = H- The four-parameter (a, 6i, 62, ^3) fits to MC data 
(empty circles) within L G [14; 128] are shown in Fig. |5|. The fit with our exponents 
(upper solid line) is relatively better at larger sizes. However, both fits (upper solid 
and dashed lines) look, in fact, quite similar, so that we cannot make unambiguous 
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conclusions herefrom. 

We have shown in Fig. ^ also the three-parameter fits to the susceptibility data 
of the 3D Heisenberg model at f3 = f3c — 0.692955 (Sec. extracted from Fig. 6 
in by a suitable linear interpolation. The MC data within L G [12; 48] are shown 
by solid circles. The fit with our exponents (r/ = 0.1, uji = 0.6/) is depicted by 
lower solid line, whereas that with the conventional RG exponents {rj = 0.0355, 
(J = 0.782 ||2l|] ) - by lower tiny-dashed line. Like in the case of 3D Ising model, 
the fit with our exponents looks slightly better, although the MC data are too 
inaccurate to make serious conclusions herefrom. Comparing the amplitudes given 
in the caption of Fig. |5|, we see that corrections to scaling are remarkably weaker in 
the 3D Heisenberg model as compared to the 3D Ising model, so that the neglected 
third-order correction in the case of the Heisenberg model could be small enough 
even at sizes somewhat below L = 12. In this aspect, it is interesting to mention 
that the fit with our exponents rj = 0.1 and uji = 0.61 (and not the other one) 
qualitatively correctly reproduces the shape of the actual In (x/i^) plot at L < 12 
where it curves upwards to meet the condition x{L = 1) = (o"^) = 1. 

For comparison, we have shown in Fig. ^ also the MC data for 3D XY model |22| 
in which, however, only x and y components of the 3-component order parameter in- 
teract with each other. One believes |22] that this model belongs to the universality 
class of the standard XY model with the number of components n = 2. Unfortu- 
nately, we have not found in the recent literature more accurate explicit data for 
n = 2 case. As we see, the actual MC data (empty boxes) at (3c evaluated approxi- 
mately (5c — 0.6444 |22] are rather scattered and, therefore, unsuitable for a refined 
analysis. Nevertheless, this is a typical situation where authors of such data make 
a very "accurate" and "convincing" estimation ^/u = 1.9696(37) or r/ = 0.0304(37) 
making a simple linear fit. However, the refined analysis given above has shown that 
even in the case of 3D Ising model, where the data are incompatibly more accurate, 
it is not so easy to distinguish between r/ = 0.0358 and r/ = 1/8. Moreover, a refined 
analysis prefer the second value which is much larger than those usually provided 
by linear fits at typical system syzes L < 48. This is particularly well seen in Fig. ^, 
where the effective critical exponent r]eff{L) of the 3D Ising model, estimated via 
the linear fit within [L; 2L], is depicted by solid circles. As we see, rjf>ff(L) tends to 
increase well above the conventional value 0.0358 (horizontal dot-dot-dashed line). 
The shape of the rjcff{L) plot is satisfactory well reproduced by a third-order poli- 
nomial in the actual scale of L~^/^. Three such kind of least-squares approximations 
(at Lrain = 9, 10, 12) are shown in Fig. ^. These fits do not provide very accurate 
and stable asymptotic values of r/. Nevertheless, they are more or less in agreement 
with our theoretical prediction 77 = 1/8 (horizontal dashed line). Besides, the values 
of rjcff are affected by the error in (ic (about 10~^) only slightly, i. e., by an amount 
not exceeding 0.001. 



6 A test for 3D Ising model with "improved" action 

Here we discuss some estimations of the critical exponents from the susceptibility 
data of 3D Ising model, reported in [^], with the so called "improved" action (i. e., 
H/T). One of the problems with the standard 3D Ising model is that corrections to 
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Figure 6: The effective critical exponent rieff{L) (solid circles) obtained by fitting the 
susceptibility data of 3D Ising model at criticality {f3 = 0.2216545) |jl^ within the inter- 
val [L;2L]. The least-squares approximations obtained by fitting the rjeff{L) data within 
[Lmin]^^ to a third-order polinomial in L~^/^ are shown by dashed {Lmm — 9), solid 
{Lmin = 10), and tiny-dashed (Lmm — 12) lines. The asymptotic value ij = 1/8 of the GFD 
theory is indicated by a horizontal dashed line. The dot-dot-dashed line represents the r] 
value 0.0358 proposed in |13|. 



scaling are strong. It has been proposed in |16| to solve this problem by considering 
a modified (spin-1) Ising model with the Hamiltonian 



H/T = -Pj2a,aj + Dj2a^ 



(11) 



where the spin takes the values 0, ±1, with two coupling constants /3 and D ad- 
justed in such a way that the leading correction to finite-size scaling vanishes for all 
relevant physical quantities (magnetization cumulant, energy per site, susceptibility, 
etc.) and their derivatives. Moreover, according to the claims in []l6| (see the con- 
clusions in ||l^), the ratios of the leading and subleading corrections are universal, 
so that not only the leading but all (!) corrections should vanish simultaneously. 

We have checked the correctness of these claims as described below. We have 
fitted the corresponding to (0) expression for In x to the susceptibility data of the 
"improved" 3D Ising model (y) with {P,D) = (0.383245,0.624235) (this is an ap- 
proximation of the critical point) given in ||l^ (Tab. 26). By fixing the exponents, 
the least-squares fit within L G [Lmin] 56] (here L = 56 is the maximal size available 
in Tab. 26 of fl^), including the leading and the subleading correction to scaling, 
provides the effective amplitudes a, 6i, and 62 depending on Lmin- We have made 
a test with the critical exponents t] = 0.0358(4), lo = 0.845(10), and u = 0.6296(3) 
proposed in [^]. These values are close to those of the usual RG expansions ||2l| 
but, as claimed in |l^, they are more accurate. According to 
expansion contains corrections like L^""^ and L^^", where n - 



the asymptotic 
1, 2, 3, . . . Thus we 
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Figure 7: The effective amplitudes IO61 (circles) and &2 (rhombs) in (||) estimated at 
fixed exponents 77 = 0.0358, uji = 0.845, uj2 = 2lji, and i/ = 0.6296 by fitting the MC 
data within L G [imi„;56]. Filled symbols correspond to 5{t,L) — 0, empty symbols - to 
S{t, L) = lO^^L^^'^. The effective amplitudes bi and &2 estimated with the critical exponents 
of our GFD theory (77 = 1/8, = 1/2) at S{t, L) = are shown by "x" and "+" , respectively. 
Lines represent the least-squares approximations by a fourth-order polinomial in L. 



have LOi = Lo and U2 = 2u!. The resulting amplitudes IO61 (Lmm) and b2{Lmin) stre 
shown in Fig. by circles and rhombs, respectively. We have depicted by filled sym- 
bols the results of the fitting with 5{t, L) = 0, assuming that the critical coupling 
f3c = 0.383245 has been estimated in [|l3| with a high enough (6 digit) accuracy. The 
data points quite well fit smooth (tiny dashed) lines within Lmin € [4; 20] , which 
means that the statistical errors are reasonably small. If the exponents used in the 
fit are correct and corrections to scaling are small indeed, then the convergence of 
the effective amplitudes to some small values is expected with increasing of Lmin- 
However, as we see from Fig. 0, the effective amplitudes tend to increase in mag- 
nitude acceleratedly as Lmin exceeds 14. A small inaccuracy in f5c value can be 
compensated by the term 6{t,L) ~ c*L^/'^ in (^), where c* = ct (cf. Eq. @)). The 
results of fitting with c* = 10~^ are shown in Fig. |^ by empty symbols. As we see, 
the expected inaccuracy in j3c of order 10~^ does not change the qualitative picture. 
The increase of the effective amplitudes indicates that either the exponents are false, 
or the asymptotic amplitudes are not small (or both). This is our argument that the 
claims in [16| about very accurate critical exponents, extracted from the 3D Ising 
model with "improved" action, are incorrect. 

For comparison, we have shown in Fig. ^ also the effective amplitudes ^i(-^^mm) 
and b2{Lmin) (by "x" and "+", respectively) estimated with the critical exponents 
of our GFD theory g 0] (r/ = 1/8, uJi = 1/2), assuming 5{t,L) = 0. The effective 
amplitudes converge to some values with increasing of L^in- These, however, are not 
the true asymptotic values, since the maximal size of the system has been eliminated 
to L = 56. 
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Figure 8: The effective amplitudes in Eq. (^) 100(lna(L) — 1) (triangles), 5bi{L) (circles), 
&2(i) (squares), and 63 (i) (rhombs) evaluated by fitting the susceptibility data of 3D Ising 
model at criticality within the interval of sizes [L; 8L] with the critical exponents rj = 1/8 
and LJi — 1/2 of the GFD theory. 

7 A test for the standard 3D Ising model 

A test with the effective amplitudes, as in Sec. ^, appears to be more sensitive tool 
as compared to the fits discussed in Sec. |5[ Since more data points are available for 
the standard Ising model, we can make even better test than that in Sec. |^. We have 
fitted all data points in Tab. 25 of |jTH| within the interval [L; 8L] to the theoretical 
expression for Inx (consistent with (g)) to evaluate the effective amplitudes a and 
hi with I = 1, 2, 3 depending on L. Exceptionally in the case if all the involved expo- 
nents are correct (exact) each effective amplitude can converge to a certain nonzero 
asymptotic value at L — > 00. In other words, if one tries to compensate the inconsis- 
tency in the exponent by choosing appropriate amplitude, then the amplitude tends 
either to zero or infinity at L ^ 00. 

We have shown in Fig. ^ the effective amplitudes In a(L) and hi{L) in the case 
of our critical exponents r/ = 1/8 and toi = 1/2. As we expected, the effective 
amplitudes converge to some nonzero values with increasing of L. This is a good 
numerical evidence that our critical exponets are true. The case with the exponents 



of ig 7/ = 0.0358(4), oji = 0.845(10), uj2 = 2uJi, and (J3 = 2 is illustrated in Fig. |. 
As we expected, the effective amplitudes of our four-parameter fit (solid symbols) 
tend to diverge with increasing of L, which shows that this set of critical exponents 
is false. One could object that, probably, the instability of the effective amplitudes 
is due to small errors in MC data. However, the amplitudes bi{L) and b2{L) of the 
more stable three-parameter fit (/ = 1,2 in (P)) behave in a similar way (see empty 
symbols in Fig. ^). Moreover, the amplitude bi{L) of the two-parameter fit, shown 
by crosses, increases almost linearly at large enough L instead of the expected (in 
a case of correct exponents) saturation like bi{L) ~ 61 + const ■ . As regards 
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Figure 9: The effective amplitudes in Eq. (^) evaluated by fitting the susceptibility data of 
3D Ising model at criticality within the interval of sizes [L; 8L] with the critical exponents 
rj — 0.0358, uji — 0.845, uj2 = 2wi, and W3 = 2 proposed in |T^. Solid symbols show the 
four-parameter fit: 506i(L) (circles), b2{L) (squares), and h'i{L) (rhombs); empty symbols 
show the three-parameter fit: 100&i(L) (circles) and 27b2{L) (squares); crosses represent 
the amplitude of the two-parameter fit, i. e., quantity 190 {bi{L) + 0.34). 



the convergence in Fig. ^ of the effective amplitudes at L — > 00, it is possible only 
if both conditions are fulfilled, i. e., the exponents are correct and the MC data are 
accurate enough to ensure stable results. Thus, in any case, the analysis in Fig. ^ 
provides rather convincing evidence that our exponents are the true ones, which 
by itself rules out the possibility that those proposed in |13| could be correct. The 
results in Figs. ^ and ^ are affected insignificantly by a small inaccuracy of about 
10~^ in the estimated /3c value. 



8 Some remarks about other numerical results 

There exists a large number of numerical results in the published literature not 
discussed here and in our previous papers 1^, . A detailed review of these results is 
given in The cited there papers report results which disagree with the values of 
the critical exponents we have proposed. However, as regards the pure Monte Carlo 
study, we are quite confident that, just like in the actually discussed case of 3D 
Ising model, the increase of system syzes and/or use of higher-level approximations 
will lead to the conclusion that fits with our exponents are better than those with 
the conventional (RG) exponents. Particularly, a careful analysis of the effective 
exponents made in Sec. |^ and ^, as well as in Sec. 6 of Q already has shown that 
the effective exponents deviate from the values predicted by the perturbative RG 
theory and converge more or less to those of the GFD theory at L — > 00. Together 
with the analysis of the experiment with superfluid '^He [^, we have presented totally 
6 independent evidences of such a behavior. 
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Formally, the finite-size effects on the obtained values of the critical exponents 
has been taken into account in many of cases considered in literature. However, 
the estimated effect strongly depends on that which kind of corrections to scaling 
is expected and included in the analysis. All the existing analysis (not counting 
our works), of course, are based on the assumption that the critical behavior of 
all physical quantities is characterised by the same correction exponent uj which is 
about 0.8 for 0(n)-symmetric models with n = 1,2, 3, 4. However, it is evident from 
the behavior of the partition function zeros of 3D Ising model that lj cannot have 
so large value. Namely, the value of (l/i^) + to should be about 2 or even smaller, 
otherwise we arrive to a rather strong and obvious contradiction with the MC data 
for the real part of the partition function zeros [Q]. The current analysis in Sec. ^ 
provides 3 independent evidences (for n = 1,2,3) that the correction exponent for 
the magnetization cumulant is remarkably smaller than 0.8. The numerical analysis 
often suggests that to ^ 1. This fact is perfectly explained by our theoretical concept: 
in some cases the amplitude of the leading correction term can be small as compared 
to that of the subleading term providing the effective correction exponent just about 
1. The value of to is crucial for an accurate correction-to-scaling analysis. If, e. g., 
we would assume that to = 0.845, then we could not arrive to a conclusion that 
r/ = 1/8 is a better choice for the 3D Ising model than rj = 0.0358(4), since all 
fits with rj = 1/8 and uj = 0.845 look relatively bad. This explains the fact that 
the usual estimations do not give r] ~ 1/8, while this is just the correct value. We 
suppose that similar problems could arise also in other cases, particularly, if one 
uses some expression for the correlation length in finite system (like in fl^), as it 
has been discussed in 0]. 

We should not forget also about purely subjective factor that any signals about 
essential inconsistency between MC data and RG predictions usually are suppressed, 
i. e., they do not apper in the published literature. There are no doubts that such 
signals exist which can be mentioned even very easily, e. g., the behavior of the 
effective critical exponents u>eff discussed in Sec. ^, or those evidences in Q which 
appear as a result of unsophisticated analysis of MC data. As a result of an un- 
critical acceptance of anything which claims to confirm with a great accuracy the 
conventional (RG) values of the critical exponents and rejection of any contraar- 
guments, the objective picture is distorted. This is the reason why almost all the 
published and reviewed papers usually claim to confirm with an almost unbelievable 
accuracy the predictions of the perturbative RG theory. It is impossible to check in 
detail all these papers, but our critical analysis in |^], and here indicates that 
many of them are, at least, inobjective. 

There exists some background for the conventional claims in the published liter- 
ature that all the usual methods give consistent results which appear to be in a good 
agreement with the predictions of the perturbative RG theory. The perturbation 
expansions of the RG theory, as well as the techniques of high- and low-temperature 
series expansion are merely not rigorous extrapolation schemes which work not too 
close to criticality. As a result, these methods produce some pseudo or effective 
critical exponents which, however, often provide a good approximation just for the 
range of temperatures not too close to Tc (critical temperature) where these meth- 
ods make sense and, therefore, agree with each other. According to the finite-size 
scaling theory, tL^f^ is a relevant scaling argument, so that not too small values of 
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the reduced temperature t are related to not too large sizes L ~ t~'^ . Therefore, one 
can understand that the MC results for finite systems often can be well matched 
to the conventional critical exponents proposed by high temperature (HT) and RG 
expansions. If, however, the level of MC analysis (i. e., the level of approximations 
used) is increased, then it turns out that the "conventional" critical exponents are 
not valid anymore, as it has been demonstrated in the current paper and in [^. It 
is because the "conventional" exponents are not the asymptotic exponents. Correct 
values of the asymptotic exponents have been found in considering suitable the- 
oretical limits instead of formal expansions in terms of In k (at criticality, where k is 
the wave vector magnitude) or In t (approaching criticality) which are meaningless at 
k and t — > 0. These formal expansions lie in the basis of the RG expansions for 
the critical exponents. The founders and defenders of the perturbative RG theory, of 
course, will try to doubt our statement that the perturbative RG method is invalid 
at criticality. But it is impossible to doubt a mathematical proof. It has been proven 
in Q that the assumption that the e-expansion works and provides correct results 
at ^ leads to an obvious contradiction in mathematics (cf. Sec. 2 in [^). This 
fact alone cannot be compensated even by an infinite number of numerical evidences 
supporting the "conventional" critical exponents coming from the RG expansions. 

Our argument, based on the current numerical analysis, is the following. We have 
proposed here a very sensitive method (i. e., a study of effective amplitudes) which 
allows to test the consistency of a given set of critical exponents with the MC data 
including several (in our case up to 3) corrections to scaling. We have applied this 
method to one of the recent and most accurate numerical data for the susceptibility 
in 3D Ising model, and have got a confirmation that our critical exponents are true. 
It would be not correct to doubt our results based on less sensitive methods and 
lower-level approximations. 

We prefer to rely just on the data of pure MC simulations becose of the following 
reasons. The so called Monte Carlo RG (MCRG) method is not free of assumptions 
related to approximate renormalization. We would like only to mention that the 
MCRG study in 1 2(: ] of 3D Ising systems of the largest (to our knowledge) available 
in literature sizes, i. e. up to L = 256, has not revealed an excellent agreement with 
the usual predictions of the perturbative RG. In particular, an estimate uj ~ 0.7 has 
been obtained [^6| which is smaller than the usual (perturbative) RG value ~ 0.8, 
but still is larger than the exact value 0.5 predicted by the GFD theory. The high- 
temperature series cannot give more precise results than those extracted from the 
recent most accurate MC data, including the actual data of |jl^, since these series 
diverge approaching the critical point. One approximates the divergent series by a 
ratio of two divergent series (Fade approximation) , but it is never proven that such 
a method converges to the exact result. It is interesting to compare the MC and HT 
estimates of the critical point for the standard 3D Ising model, i. e., /3c — 0.2216545 
(MC) fill and /3c = 0.221659 + 0.000002/ - 0.000005 (HT) [||]. It is clear that the 
MC value is more accurate: if we look in ||l^, where the estimation procedure is 
well illustrated, we can see that (3c is definitely smaller than 0.221659, and the error 
seems to be much smaller than the difference between both estimates 0.0000045. 
As we have mentioned already, our independent tests suggest that the error of the 
actual MC value is about 10~^. 
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9 Conclusions 



In summary of the present work, we conclude the following. 

1. The leading correction-to-scaling exponent iv in 0(n)-symmetric models on a 
three-dimensional lattice has been tested by analysing the recent Monte Carlo 
(MC) data. These tests have shown the incorrectness of some claims that uj 
has a very accurate value 0.845(10) at n = 1. A selfconsistent infinite volume 
extrapolation yields row estimates uj 0.547, uj 0.573, and uj « 0.625 at 
n = 1, 2, and 3, respectively, in approximate agreement with the corresponding 
exact values 1/2, 5/9, and 3/5 predicted by our recently developed GFD theory. 

2. Considering the susceptibility data for 3D Ising model at criticality (Sec. |5|), we 
conclude that the fits made with our (GFD) critical exponents systematically 



improve relative to those made with the exponents given in |13|, as the system 
sizes grow and the approximation order increases. 

3. The numerical analysis of the effective critical exponents in Sec. I and |, as 
well as in Sec. 6 of Q has shown that the effective critical exponents deviate 
from the values predicted by the perturbative RG theory and converge towards 
those of the GFD theory at L ^ oo. The same behavior has been observed 
in the experiment with superfluid ^He discussed in Q. Totally, these are 6 
independent evidences of such a behavior, suggesting that the above examples 
are not occasional or exceptional, but reflect a general rule. 

4. Different sets of critical exponents (one provided by GFD theory, another pro- 
posed in |]l3|) predicted for the 3D Ising model have been tested by analysing 
the effective amplitudes (Sec. ^ and ^. While the usual fits of the suscepti- 
bility data do not allow to show convincingly which of the discussed here sets 
of the critical exponents is better, this method clearly demonstrates that the 



conventional critical exponents rj = 0.0358(4) and uj = 0.845(10) [13| are false, 
whereas our (GFD) values rj = 1/8 and a; = 1/2 are true. 
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